Structure of AlSb(OOl) and GaSb(OOl) Surfaces Under Extreme Sb-rich Conditions 
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We use density-functional theory to study the structure of AlSb(OOl) and GaSb(OOl) surfaces. 
Based on a variety of reconstruction models, we construct surface stability diagrams for AlSb and 
GaSb under different growth conditions. For AlSb(OOl), the predictions are in excellent agreement 
with experimentally observed reconstructions. For GaSb(OOl), we show that previously proposed 
model accounts for the experimentally observed reconstructions under Ga-rich growth conditions, 
but fails to explain the experimental observations under Sb-rich conditions. We propose a new 
model that has a substantially lower surface energy than all (nx5)-like reconstructions proposed 
previously and that, in addition, leads to a simulated STM image in better agreement with ex- 
periment than existing models. However, this new model has higher surface energy than some of 
(4x3)-like reconstructions, models with periodicity that has not been observed. Hence we conclude 
that the experimentally observed (1x5) and (2x5) structures on GaSb(OOl) are kinetically limited 
rather than at the ground state. 

PACS numbers: 68.35.Bs, 68.35.Md, 68.37.Ef, 68.47.Fg 
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I. INTRODUCTION 

The surfaces and interfaces of III-V semiconductors 
constitute some of the most important components of the 
semiconductor industry. For example, III-V heterostruc- 
ture quantum wells are key components in a wide range of 
optical and high-frequency electronic devices, including 
field-effect transistors^, resonant tunneling structures^, 
infrared lasers'^, and infrared detectors^. Many of these 
devices require extremely sharp and clean interfaces. For 
this reason, an understanding of the atomic-scale mor- 
phology of III-V semiconductor surfaces is critical to con- 
trolling the growth and formation of their interfaces^i^. 

It is generally accepted that the surfaces of III-V semi- 
conductors should reconstruct in such a way that the 
number of electrons is exactly enough to doubly occupy 
all orbitals on electronegative (V) atoms, leaving all or- 
bitals on electropositive (III) atoms unoccupied. This 
guiding principle, known as the electron-counting model 
(ECM), has been used to screen candidate structural 
models of many observed reconstructions on the surfaces 
of III-V semiconductors^i^'^iiSiiii. In practice, however, 
not all experimentally realized reconstructions follow this 
principle. For example, under Sb-rich growth conditions, 
GaSb(OOl) forms surface reconstructions that are weakly 
metallic and hence violate the ECM^^, even though the 
closely related AlSb(OOl) surface forms insulating recon- 
structions that satisfy i1>i^. The nature of reconstruc- 
tions that violate the ECM, and the underlying reasons 
for their stability, are thus important for understanding 
III-V surfaces in general. 

In this article we explore theoretically a large num- 
ber of judiciously chosen candidate reconstructions on 
GaSb(OOl) and AlSb(OOl). We find that as the growth 



conditions are varied between Sb-poor and Sb-rich, 
the predicted sequence of stable reconstructions for 
GaSb(OOl) is exactly analogous to those of AlSb(OOl). 
Experimentally, however, the picture is more compli- 
cated. In the Sb-poor limit, the observed GaSb(OOl) re- 
construction is indeed analogous to that of AlSb(OOl). 
On the other hand, in the Sb-rich limit, the exper- 
imentally observed reconstructions for GaSb(OOl) and 
AlSb(OOl) are different. Moreover, in this limit the pre- 
dicted and observed reconstructions are in good agree- 
ment only for AlSb(OOl), while for GaSb(OOl) there re- 
mains an unresolved discrepancy between theory and ex- 
periment. 

Experimentally, the Sb-terminated AlSb(OOl) surface 
evolves through the sequence Q!(4x3) — > /3(4x3) — * 
7(4x3) c(4x4) as the growth condition is changed 
from low Sb flux (or high substrate temperature) to high 
Sb flux (or low temperature)^. All of these reconstruc- 
tions are insulating, and are well accounted for by struc- 
tural models proposed in the literature that satisfy the 
ECM. 

Of particular interest here is the Sb-rich AlSb(OOl)- 
c(4x4) reconstruction, analogous to the As-rich 
GaAs(001)-c(4x4) reconstruction, which is observed on 
AlSb but not GaSb. In contrast to both AlSb and 
GaAs, the GaSb(OOl) surface does not exhibit a sta- 
ble, insulating c(4x4) reconstruction under similar — or 
any other — growth conditions. Instead, it forms (nx5)- 
like reconstructionsi^iiiii^iiSiiil. Structural models pro- 
posed in the literature for these (nx5)-like reconstruc- 
tions violate, by design, the ECM and consequently are 
weakly metallic^^. Simulated scanning tunneling mi- 
croscopy (STM) images based on (2x10) and c(2xl0) 
models closely resemble the experimental imagesi?. As 
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a result, these models have been generally accepted as 
describing the GaSb(OOl) surface under Sb-rich growth 
conditions. Nevertheless, we show below on energetic 
grounds that these models are unlikely to be correct. 
Specifically, we find their calculated surface energy to 
be significantly higher than GaSb(001)-c(4x4) for any 
plausible value of Sb chemical potential. However, since 
the experimentally observed reconstruction of GaSb(OOl) 
does not have c(4x4) periodicity, this model cannot be 
correct either. Thus a definitive structural model remains 
to be found. 



II. METHODS 

The basic structural models we considered are taken 
from the literature and are shown in Figs. [T] and [2l Sur- 
faces that satisfy the ECM are generally semiconducting, 
while those that do not may be metallic. The degree to 
which a given surface satisfies the ECM can be measured 
by the excess electron count, /S.v, which we define here as 
the difference between the number of available electrons 
and the number required to satisfy the ECM, per (1x1) 
surface unit cell. Excess electron counts for the structural 
models in Figs. [T] and [2] are tabulated in Table HI 

To compare the surface energies of reconstruction mod- 
els with different periodicities and stoichiometrics, we 
consider the surface energy per unit area, 

7 = Esurf/A = {Etot ~ miiMni - '^vMv)Mi (1) 

where i?tot is the total energy of a reconstructed sur- 
face, of area A, containing nm group-Ill and rtv group- 
V adatoms in excess with respect to the bulk-truncated, 
Sb-terminated AlSb(OOl) or GaSb(OOl). The atomic 
chemical potentials /i' are more conveniently expressed 
in terms of excess chemical potentials /i, relative to the 
energy per atom in the ground-state elemental phases: 
/i' = /i^""^ + ji. Assuming the surface to be in thermo- 
dynamic equilibrium with the bulk, the III and V chem- 
ical potentials are related by -|- /iy — ^Hf, where 
AHf = fi^^\ - (/i^fji + /i^""^) is the formation enthalpy 
of the bulk III-V crystalline phase^^- (note that AHf is 
intrinsically negative). Eq. ^ can then be rewritten to 
show more clearly the dependence of 7 on the surface 
stoichiometry and chemical potential: 

7 = 70 + AivAe. (2) 

Here 70 = {Et - S.^b) - Mm-v0in + My A© is indepen- 
dent of the chemical potentials, and A8 = Om — Qy — 
(niii — ny)/^ is the deviation of the surface stoichiome- 
try from its bulk value. The dependence of 7 on chemi- 
cal potential is given entirely by the second term. Note 
that iiy is intrinsically negative, and can take values in 
the range AHf < fiy < 0. Hence, Eq. ([2]) reflects in a 
simple way that Ill-rich reconstructions {Qm > Qy) are 
favored under Ill-rich conditions (/iy — > AHf), V-rich re- 
constructions (0v > 0iii) are favored under V-rich con- 




Y(4x3) c(4x4) 



FIG. 1: (color online) Reconstruction models proposed for the 
AlSb(OOl) or GaSb(OOl) surfaces with (4x3) and (4x4) peri- 
odicities. The first two upper layers are shown in a top view. 
Smaller white circles represent Sb atoms in the top layer of 
the underlying Sb-terminated AlSb(OOl) or GaSb(OOl) sur- 
face. Larger circles represent AI or Ga (black) and Sb (white) 
adatoms. The unit cells are shown in light blue. 



ditions (/iy — s- 0), and for stoichiometric reconstructions 
(Oy = Qui) 7 does not depend on chemical potential. 

To compute the total-energy contribution, 70, to the 
surface energy we performed first-principles calculations 
using density-functional theory (DFT). The calculations 
were performed within the local-density approximation 
(LDA)^~!^ using ultrasoft pseudopotentials^^'^^'^'^. We 
used a standard supercell technique, modehng the (001) 
surface with a slab consisting of four bilayers and 10 A 
of vacuum. Atoms in the bottom bilayer were fixed at 
their bulk positions, while all other atoms are allowed 
to relax until the rms force was less than 0.005 eV/A. 
The bottom layer (either Ga or Al) was passivated with 
pseudohydrogens. A plane-wave cutoff of 300 eV was 
used, and reciprocal space was sampled with a density 
equivalent to at least 192 fc-points in the (1x1) surface 
Brillouin zone. To define the III-V formation enthalpy 
AHf from the bulk chemical potentials /i*^""*, separate 
DFT calculations were performed for the elements in 
their ground-state phases: Ga in the a-Ga structure, Al 
in the face-centered cubic structure, Sb in the rhombo- 
hedral structure, and both AlSb and GaSb in the zinc 
blende structure. 
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TABLE L Electron count for different reconstructions of GaSb(OOl) surface. Tlie excess electron count per (1x1) surface unit 
cell is defined as Az/ — {h — m)/A where n is the number of available electrons and m is the number of required electrons to 
satisfy the ECM in the excess of Sb-terminated GaSb(OOl). A is the area of the surface unit cell in terms of the (1x1) surface 
unit cell, rii is number of adatoms of species i in excess with respect to the Sb-terminated GaSb(OOl) and Qi — rii/A is the 
coverage of adatoms of species i. The relative 7 values, in eV per (1x1) surface unit cell, are given with respect to that of 
a(4x3). 
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represent the second layer Sb adatoms. 

FIG. 3: (color online) Reconstruction models with a single 
substitution of Sb atoms by Ga atoms. See Fig. [T] for color 
III. RESULTS AND DISCUSSIONS schemes. 



The resulting relative surface energies for AlSb(OOl) 
and GaSb(OOl) are shown in Figs-EJa) andElJb), respec- 
tively, for the eight models considered here. For each 
model, the surface energy is linear in /iy, with the slope 



given by AO. 

For AlSb(OOl) the predicted stable reconstructions, 
and their energetic ordering, are in excellent agreement 
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FIG. 4: (color online) Reconstruction models with a double substitution of Sb atoms by Ga atoms. See Fig.[T]for color schemes. 



with experiment. Proceeding from Sb-poor to Sb-rich 
conditions, the predicted sequence is q;(4x3) — > /3(4x3) 
7(4x3) c(4x4), as reported previously This is 
the same sequence observed experimentally^"^. Moreover, 
7(4x3) is predicted to exist only over a very narrow range 
of /isb, in agreement with experiments^. 

For GaSb(OOl) the predicted sequence is qualitatively 
the same as for AlSb(OOl), although the c(4x4) is only 
predicted to be stable for values of /isb above the ther- 
modynamically allowed limit of zero. Experimentally, 
however, the situation is quite different. As reported 
previously, neither the 7(4x3) nor the c(4x4) phase is 
observed for any growth conditior^i^. Instead, under Sb- 
rich conditions, only the (1x5) and (2x5) periodicities 
have been observedi^. 



Righi et al. suggested /i0(4x3), shown in Fig.[Tl as the 
model for GaSb(OOl) surface under these conditions^!. 
Our calculation indeed shows that it is energetically as 
favorable as 7(4x3), as shown in Fig.[5l^b) and in TableU 
However, ft,0(4x3) must be rejected as it has a wrong 
periodicity. 

In order to explain the experimentally observed (1x5) 
and (2x5) structures on GaSb(OOl) surface, we studied a 
large number of structures based on variations of c(2 x 10) 
and (2x10). We note that c(2xl0) violates the ECM 
substantially {Av = 0.3) and substitution of Sb atoms in 
the top layer of the underlying Sb-terminated GaSb(OOl) 
surface by Ga atoms can lower the excess electron count. 
Fig. [3] shows the possible reconstructions when a single 
Sb atom is replaced by a Ga atom. We use the naming 
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FIG. 5: (color online) (a) Surface stability phase diagram for 
AlSb(OOl) surface. The relative surface energy [Eq. ([2])] is 
plotted as a function of the Sb chemical potential relative to 
its corresponding bulk value. Dotted vertical lines mark the 
thermodynamically allowed range of fisb- ^Hf is the heat of 
formation for AlSb. (b) Surface stability phase diagram for 
GaSb(OOl) surface; AHf is the heat of formation for GaSb. 



convention of six to denote a "single substitution". As 
shown in TablelH all six reconstruction indeed have lower 
excess electron counts. 

For completeness, we also considered reconstructions 
resulting from double substitution of Sb atoms by Ga 
atoms as shown in Fig. [4] More substitutions, however, 
were not found to be energetically favorable: Table H] 
shows that the surface energy of these structures are 
higher than that of six reconstructions. We note that 
for these double substitutions the excess electron counts 
Av are negative, indicating a deficit of electrons relative 
to the ECM. 

One of the most energetically favorable structures 
having the correct periodicity is sla-c(2xl0), shown 
in Fig. [3] sla-c(2xl0) has two clear advantages over 
c(2xl0). First, the surface energy for sla-c(2xl0) is 



FIG. 6: Filled-state STM images of GaSb(OOl) with (1x5) pe- 
riodicity, (a) Experimental STM image; a c(2 x 10) unit cell 
is indicated, (b) Simulated STM image of sla-c(2xl0). This 
image shows the asymmetries in the intensities of the cur- 
rent density from two atoms of the horizontal dimers, which 
was not captured in the simulated STM image of c(2xl0). 
Compare with Fig. 3(d) of Ref. [H. 



lower than that of c(2xl0) by more than 50 meV per 
(1x1) unit cell. Second, as shown in Fig. [SI the simu- 
lated STM image for sla-c(2xl0) is in a better agree- 
ment with the experiment image, in that it reproduces 
the left-right asymmetry within the surface Sb dimersi^. 
Furthermore, as shown in Table IH this model violates 
the ECM and thus is predicted to be weakly metallic, 
as observed in tunneling spectroscopy^. Therefore, the 
previously proposed model c(2xl0) is unlikely to be the 
experimentally realized structure. 

However, the calculated surface energy of sla-c(2xl0) 
is higher than that of 7(4x3), as shown in Tableland 
Fig. IHb). Likewise, (2x10), the structural model gen- 
erally accepted for the surface with (2x5) periodicity, 
is the least energetically favorable structure among the 
eight structures of Table [J On the other hand, 7(4x3), 
the most energetically favorable structure among all the 
structures considered in this study, has a periodicity that 
has not been observed experimentally to date. These 
facts leave us with two possible conclusions: either the 
correct structural model remains undiscovered, or the 
experimentally obtained surface is not the ground-state 
structure. 

The latter possibility, a kinetically limited surface, 
bears closer consideration. For example, there may be an 
activation barrier to forming the mixed dimers on GaSb 
that cannot be overcome within the growth temperatures 
and times used here. Indeed, to stabilize these surfaces 
during the growth, one must go from active growth with 
both Ga and Sb flux at ~500 °C, to room temperature 
and no flux while trying to stabilize the surface. This 
process typically involved simultaneously lowering the 
temperature while turning off the Ga and then lower- 
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ing the Sb flux. The surface cannot be annealed, be- 
cause that would drive off Sb and create (nx3) recon- 
structions. These considerations lead us to propose that 
the sla-c(2xl0) structure is the most likely model for 
the observed GaSb(OOl) surface as created under Sb-rich 
growth conditions and subsequently stabilized under vac- 
uum. 



IV. SUMMARY AND CONCLUSIONS 

We have performed ab initio calculations on the surface 
energy and atomic structure of AlSb(OOl) and GaSb(OOl) 
surfaces with various reconstructions. Surface stability 
diagrams for a large number of reconstruction models 
are constructed under different growth conditions. For 
AlSb(OOl), we confirmed that the predictions of the cur- 
rently accepted models are in good agreement with ex- 
perimentally observed reconstructions. For GaSb(OOl), 
we showed that previously proposed model accounts for 
the experimentally observed reconstructions under Ga- 
rich growth conditions, but fails to explain the experi- 
mental observations under Sb-rich conditions. Therefore, 



we propose sla-c(2xl0) as a better alternative to exist- 
ing models for GaSb(OOl) under extreme Sb-rich growth 
conditions. Our calculations show that sla-c(2xl0) has 
a substantially lower surface energy than all (nx5)-like 
reconstructions proposed previously and, in addition, it 
leads to a simulated STM image in better agreement with 
experiment than existing models. However, sla-c(2xl0) 
has higher surface energy than 7(4x3), a model with pe- 
riodicity that has not been observed. Hence we conclude 
that the experimentally observed (1x5) and (2x5) struc- 
tures on GaSb(OOl) are not the ground-state structure, 
but kinetically limited ones. 
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